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ABSTRACT 

We analyze radial velocities for a sample of 31 Class I and flat spectrum protostars in Taurus- 
Auriga, p Ophiuchi and Serpens for evidence of the global dynamical state of extremely young stellar 
populations buried within parental molecular clouds. Comparing the radial velocity of each protostar 
to that of the local CO gas, we are able to constrain the one dimensional radial velocity dispersion 
of Class I and flat spectrum objects to ~2.5 km sec -1 or below. This upper limit to the protostellar 
velocity dispersion is consistent with the velocity dispersions of surrounding CO gas which we measure 
to be ~ 1.4 km sec -1 , suggesting that the motions of protostars and local CO gas are dynamically 
linked and dominated by the gravitational potential of the molecular cloud. However, the upper 
limit on the protostellar velocity dispersion could still allow for slightly inflated motions of protostars 
relative to the local molecular gas. Four of the protostars analyzed appear to have velocities more than 
3 a (7.5 km sec -1 ) away from the central local CO gas velocity while showing spectroscopic indicators 
of youth and accretion such as H2 emission, HI Br 7 emission, or K band continuum veiling. These 
radial velocity outliers may represent protostellar spectroscopic binaries or ejected cluster members. 
Subject headings: infraredrstars — stars:formation — stars:low-mass — stars:pre-main-sequence — 
stars: kinematics 



1. INTRODUCTION 

The dynamical state of embedded protostellar popula- 
tions is of crucial interest for the star formation process. 
Until now, however, a diagnosis of the dynamical state 
of young (Class I and flat spectrum) protostellar popula- 
tions has been largely inaccessible to direct observational 
study. In this study we investigate the kinematics of low 
mass protostars and interpret the results within existing 
models of the star formation process. 

Theoretical investigations have recently identified dy- 
namical interactions between accreting protostars as a 
possible critical factor in the star formation process; this 
has been suggested both through analytic arguments 
l)R,einurth fc Clarke! 120011) and N-body sim ulations of 
star f ormation in a small cluster environment iBate et alJ 
2003) . These studies indicate that gravitational encoun- 
ters between protostars can separate a heavily accret- 
ing protostar from its mass reservoir, preventing further 
mass accretion and stunting its growth. These interac- 
tions would predominantly affect the smallest members 
of an association, whose lower masses make them more 
susceptible to accelerations from gravitational interac- 
tions with nearby neighbors. Accordingly, it has been 
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suggested that dynamical interactions or ejections from 
an unstable multiple system may play a large role in 
shaping the mass spectrum in the substellar regime. 

Observational clues to the importance of the ejection 
mechanism in the context of the larger star formation 
process, however, are mixed. Young brown dwarfs have 
been observed to po ssess infrared excess e s indicative of 
circu mstellar disks ( Wilki ng et al.1 Il99?l iMuench et alJ 
I2001D and signatures of activ e accretion have been ob- 
served in their spect ra as well (jJavawardhana et all2003l 
iLuhman et alJl2003|) . The predominance of steady mass 
accretion in isolated young brown dwarfs would seem 
to suggest any possible ejection mechanism must either 
not completely disrupt the accretion reservoirs of brown 
dwarfs (the very objects which should be most suscepti- 
ble to disruption) or rare enough not to produce a siz- 
able fraction of the br own dwarf population. However, 
iMuzerolle et alJ l)200 3) find that analysis of Ha emission 
from young brown dwarfs indicates an accretion rate of 
M cx M~ 2 between 1 and 0.05 Mq, and that only 3 of 13 
brown dwarfs in their sample appear to be strongly ac- 
creting. These results may indicate that isolated brown 
dwarf accretion could be weaker than we might expect 
from linearly scaling T Tauri accretion to lower masses, 
possibly as a result of prior disturbances by interactions. 

A more direct probe of the ejection paradigm would be 
to directly measure the kinematic signature of ejection in 
sites of ongoing star formation. We expect significant dy- 
namical interactions between protostars will inflate the 
intrinsic velocity dispersion of the protostars above that 
of the local molecular gas, while a star formation pro- 
cess in which dynamical effects have little influence will 
result in a protostellar velocity dispersion similar to that 
of the local molecular gas from whence the protostars 
formed. In short, an increase in the protostellar radial 
velocity dispersion above that of the local molecular gas 
serves as a proxy for the degree to which interactions 
affect the star formation process. This intuition is but- 
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tressed by the results of numerical simulations of the star 
formation process, which show stellar interactions can in- 
crease the protostellar velocity dispersion by a factor of 
three above the velocity dispersion of the m olecular cloud 
from which they formed l)Bate et al.ll200.1l) . Precise mea- 
surements, however, are required to test this model, as 
N-body simulations predict velocity dispersions even in 
such interacting environments at the 1-2 km sec -1 level. 

The dynamic state of protostellar populations is also 
interesting in the context of cluster formation and dis- 
ruption. The majority of young, forming stars in the 
Milky Way are to be found in cluster environments 
ijLada fc Ladal l2003lh However, the stellar population 
of the Milky Way galaxy is dominated not by clusters 
but by the diffuse field population making up the thin 
and thick disk populations, presu mably from stars which 
formed in now dissipated clusters l)Maiew ski 1993). Ex- 
plosive gas removal is thought to efficiently disrupt the 
majority of young clusters on timescales of 5-10 Myrs. 
Those clusters which do survive the dispersal of their 
parental gas cloud eventually fall prey to evaporation 
and tidal heating, process es which operate o n longer 
timescales of 10 8 " 9 years (|Lada fc Ladal 12003ft . These 
effects are thought to leave measurable kinematic sig- 
natures on the populations they effect. Mass loss from 
clusters, for example, is thought to result in a decreased 
velocity dispersion for the remaining bound members; 
the ratio of the initial and final velocity dispersions are 
correlat ed with the ove rall star formation efficiency in the 
cluster (Mathicu 1983). Measuring the dynamic state of 
protostellar populations can provide insight into mech- 
anisms of cluster disruption by providing constraints on 
the initial velocity dispersion of newly formed stars for 
numerical simulations of cluster disruption, as well as by 
probing cluster evolution with a comparison of protostel- 
lar populations to more evolved bound clusters. 

The velocity dispersion of protostars in embedded clus- 
ters has not yet been observationally determined. The ra- 
dial velocity dispersion has been studi ed for more evo lved 
T Tauri stars beginning in the w ork of[ Herbig] <|197 7[) an d 
continuing through the wo rk of IJones fc Herbid l)1979|) . 
lHartmann et all l)1986ft and lDubath et al.l (|1996^ . These 
studies find no difference above the errors (typically ~ 
0.5 km/sec) between the mean stellar velocity and the 
central velocity of the local gas, and measure 1-D radial 
velocity dispersions of 2 km sec -1 or less. A study of this 
type has not yet been performed for protostellar popula- 
tions, primarily because the spectroscopic resolution and 
sensitivity needed to measure radial velocities at the 2-3 
km sec -1 level has only recently become accessible in the 
infrared where these very young, highly embedded and 
reddened objects are easier to detect. 

Radial velocity measurements also serve to constrain 
the fraction of stellar systems which are members of 
close multiple systems. Recent near infrared imaging 
surveys of Class I and flat spectrum sources in the 
Taurus, p Ophiuchi, fc Serpens star for mation regions 
ijHaisch et al.ll2004t IDuchene et al.ll2004j) show compan- 
ion star fractions of ~ 20 — 25% for separations be- 
tween 300 and 2000 AU. The companion star frac- 
tion of embedded protostars is thus significantly ele- 
vated above that detected for solar type main sequence 
stars in the same range of sepa rations (12%: iMathievJ 
11994 iDuauennov fc MavoiIl 99lL as also found for more 



evolve d T Tauri stars in the same clouds. IDuchene et alJ 
( 2004) also find that Class and I sources with millime- 
ter emission from protostellar envelopes have higher com- 
panion star fractions (38%) than Class I/II/III objects 
without significant millimeter envelope emission (22%). 
These results suggest the multiplicity fraction of stars 
may decline through the pre-main sequence phase from 
high fractions in the most heavily embedded phases to 
lower fractions on the main sequence. The observa- 
tions analyzed here are sensitive to binaries with large 
orbital velocities, sampling bina ries with much sm aller 
separations than the studies of lHaisch et alJ i|2004l) and 
IDuchene et alJ 1)2004 ). useful for probing possible effects 
of orbital evolution between the protostellar, T Tauri and 
main sequence phases. 

In this study, we investigate the dynamical motions 
of Class I and flat spectrum protostars by studying the 
radial velocity dispersion of protostars about the local 
CO gas velocity. In §2 we describe the observations used 
to derive radial velocities for Class I and flat spectrum 
protostars, as well as the observations which allow us to 
determine the systemic velocity of the CO gas associ- 
ated with each star formation region. Section 3 contains 
our analysis of the resultant radial velocity dispersion 
for Class I and flat spectrum sources, while §4 constrains 
the possible source of four prominent protostellar radial 
velocity outliers. Section 5 discusses the implications of 
these results in the context of the standard picture of 
low mass star formation, and a summary of this work is 
contained in §6. 

2. OBSERVATIONS & VELOCITY DETERMINATIONS 

This study utilizes data from two distinct sources; the 
near IR spectral surv ey of Class I and flat spe ctrum pro- 
tostars presented bv iDoppmann et al] l)2005f) (hereafter 
D0 5), and the CO su rvey of the Milky Way described 
bv lDame et"aTl l|200lD . Though we refer the reader to 
the publications cited above for a full and complete dis- 
cussion of the observations conducted by each survey, 
we provide below a brief summary of those observational 
details most relevant to the work presented here. 

2.1. Protostellar Observations and Derived Radial 
Velocities 

The protostellar radial velocities used in this study 
were derived using sprectra from the Class I and flat spec- 
trum survey conducted by D05. Class I and flat spectrum 
sources were observed by D05 over the course of 12 nights 
between May 2000 and June 200 3 using the near infra red 
echelle spectrograph NIRSPEC <|McLean et al.lll998j) on 
the 10 m Keck II telescope. Spectra were acquired with a 
0.58" (4 pixel) wide slit (allowing R ~ 17,000), utilizing a 
1024x1024 pixel InSb detector array with the NIRSPEC- 
7 blocking filter, under good (~ 0.6") seeing. The NIR- 
SPEC gratings were configured such that 4 prominent 
sets of spectral lines could be simultaneously observed in 
non-continuous orders: 1) the 2.1066 /iin Mg and 2.1099 
[Mm Al lines, 2) the 2.16609 /jm HI 7 line, 3) the 2.2062 
and 2.2090 /im Na doublet (blended with lines of Si and 
Sc) and the 2.2233 ^m H 2 line, and 4) the 2.2935 /mi 
CO bandhead. We have selected for our analysis pro- 
tostars with detected absorption features located in the 
p Ophiuchus, Taurus-Auriga or Serpens star forming re- 
gions and observed by D05 in July 2001 or later, as the 
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May 2000 data presented by D05 has much larger wave- 
length calibration uncertainties. These criteria resulted 
in a sample of 31 protostars suitable for analysis in this 
work. The set of MK spectral type standards observed 
by D05 have also been included in this analysis to provide 
estimates of observational uncertainties. 

A discussion of the methods used to reduce spectral 
data is presented in full by D05. In summary, all spec- 
tra were flat-fielded, cleaned, sky subtracted, extracted, 
wavelength calibrated and corrected for telluric absorp- 
tion using standard IRAF packages and procedures. The 
radial velocity of each star was determined through cross 
correlation (using IRAFs FXCOR package) with spectral 
standards observed by D05 and with well measured radial 
velocities ( Nid ever et al.ll2002]) . The targets examined in 
this work range in spectral type from late G to mid M; 
given the onset of significant continuum structure for ob- 
jects of types MO and later, we divided the targets into 
pre and post MO sub-samples for the purposes of cross- 
correlation. Targets with spectral types earlier than MO 
were cross correlated with template spectra of HD166620 
(K2), HD201091 (K5) and HD219134 (K3). Targets with 
spectral types of MO or later were cross correlated with 
template spectra of Gliese 806 (M2), Gliese 4281 (M6.5) 
and VB8 (M7). Cross correlation was performed on 3 
echelle orders containing Mg, Al, Na and CO absorption 
features; each object was compared to all three applicable 
templates using all orders with detected absorption fea- 
tures. A mean heliocentric velocity for each object was 
then found as the average of 9 independently derived ve- 
locity estimates (from 3 orders in each of 3 standards) in 
all but two cases. Disk emission in EC38's CO order re- 
duced the number of independent velocity estimates to 6, 
while WL12 produced extremely poor velocity estimates 
via cross correlation, so we have adopted D05's radial 
velocity estimate for this object. The radial velocity of 
these protostars with respect to the local standard of rest 
(LSR) is presented in Table 1, after corrections for the 
solar motion with respect to the LSR were applied using 
IRAFs RVCORRECT task. 

To test the accuracy of the radial velocity determina- 
tion resulting from this procedure, we have compared our 
derived heliocentric radial velocities with reported values 
for 15 spectral standards observed by D05 with accu- 
rately determined radia l velocities (errors < 2 km sec" 1 ; 
[ Mohantv fc Basril 120031: iNidever et al.ll2002t iGizis et all 
l2002tlde Medeiros fc Mavorll 1999ft . Figure 1 displays the 
residuals obtained when the radial velocity reported in 
the literature is subtracted from our derived value. We 
assume that the spread of non-zero residuals in Figure 
1 is due primarily to observational error - that is, we 
believe Figure 1 represents the expected observed radial 
velocity distribution for a population of stars with iden- 
tical radial velocities (a = km/sec) and observational 
errors typical of the limits of the observations presented 
by D05. A best fit guassian to this distribution, shown 
in binned form in Figure 1 as a dashed histogram, re- 
sults in a gaussian width (a) of 1.28 km sec -1 . We thus 
conservatively adopt 1.5 km sec -1 as an estimate of the 
uncertainty in our velocity determinations for slowly ro- 
tating objects with strong absorption features. 

2.2. Molecular Cloud Observations and Derived Radial 

Velocities 



TABLE 1 

Radial Velocities of Class I and Flat Spectrum Protostars 



Target name 


V; sr 
(km/sec) 


CO V Ur 
(km/sec) 


di f f 

(km/sec) 


CO-stellar dist 
(beam widths) 


04016+2610 


-10.42 


6.83 


-17.25 


0.56 


04489+3042 


-18.27 


6.83 


-25.10 


0.61 


04295+2251 


3.83 


5.53 


-1.70 


0.48 


Haro-6-28 


3.64 


5.53 


-1.89 


0.26 


04108+2803B 


6.59 


6.83 


-0.24 


0.43 


DG-Tau 


8.20 


6.83 


1.37 


0.45 


GV-Tau S 


-6.22 


6.18 


-12.40 


0.54 


04264+2433 


6.72 


6.18 


0.54 


0.40 


04158+2805 


5.47 


7.48 


-2.01 


0.43 


04181+2655 


7.55 


6.83 


0.72 


0.39 


L1551-IRS5 


7.89 


6.83 


1.06 


0.55 


IRS43 


3.01 


3.74 


-0.73 


0.43 


WL12 


-0.61 


3.48 


-4.09 


0.19 


IRS67 


3.26 


4.78 


-1.52 


0.78 


YLW16A 


4.41 


3.74 


0.67 


0.57 


162636 


3.97 


3.22 


0.75 


0.42 


GY21 


2.65 


3.22 


-0.57 


0.83 


GY224 


2.30 


3.74 


-1.44 


0.82 


VSSG17 


7.06 


3.48 


3.58 


0.68 


WL3 


3.71 


3.48 


0.23 


0.84 


WL17 


4.88 


3.48 


1.40 


0.89 


IRS51 


6.30 


3.74 


2.56 


0.12 


GY197 


6.52 


3.48 


3.04 


0.70 


GY91 


2.91 


3.48 


-0.57 


0.85 


CRBR12 


-4.60 


3.22 


-7.82 


0.62 


EC125 


5.07 


7.91 


-2.84 


0.92 


EC38 


12.88 


8.56 


4.32 


1.22 


EC91 


5.66 


7.91 


-2.25 


0.77 


EC92 


13.02 


7.91 


5.11 


0.82 


EC94 


9.70 


7.91 


1.79 


0.79 


EC129 


10.51 


7.91 


2.60 


0.71 


EC53 


13.36 


8.56 


4.80 


1.26 



The radial velocity dispersion of the protostars in each 
individual star forming region can be calculated directly 
from the Vlsrs reported in Table 1. As we have only 
observed ~ 10 protostars in each star formation region, 
individual radial velocity dispersions for each star for- 
mation region would be very sparsely sampled, possi- 
bly returning spurious results. To improve the statistical 
leverage of our analysis, we have chosen to investigate 
the motions of protostars relative to the central velocity 
of the local molecular gas. Assuming that the dynamical 
behaviors in the Taurus- Auriga, Serpens and p Ophiuchi 
star formation regions are equivalent to first order, study- 
ing protostellar velocity residuals will allow us to combine 
results derived from protostars in different star formation 
regions, generating a snapshot of the dynamical state of 
newly formed stars. The dispersion of protostellar radial 
velocities about the local gas velocity can also be com- 
pared to the velocity dispersion of the gas itself, allowing 
investigations of the extent to which stars and gas re- 
main dynamically linked during the Class I/flat spectrum 
phase. This enhanced statistical power comes at the ex- 
pense of our ability to investigate differences in dynamics 
between star formation regions. Additionally, if ejection 
and evaporation are already important dynamical mech- 
anisms by the Class I/flat spectrum stage, we might ex- 
pect protostellar surveys to preferentially detect candi- 
dates whose motion has carried them to regions of the 
cloud with lower optical depths, as viewed from Earth. 
This effect would be detected by a protostellar radial ve- 
locity distribution with a peak at lower (bluer) velocities 
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Spectral type standard velocity residuals (Obs - Lit) 



Protostar residuals 
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Gauss center = -0.627161 
Gauss width =1.27924 
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FlG. 1. — Histogram display of residuals (solid line) when heliocentric 
velocities derived from NIRSPEC spectra of spectral type standards are 
compared to accurate (errors < 2 km sec -1 ) determinations from the 
literature for the same stars. The best fit guassian to the histogram, 
binned in the same manner, is shown as the dotted line. The gaussian 
a provides an indication of the level of accuracy with which we can 
assign velocities {a ~ 1.3 km sec -1 ). 



than the associated CO gas, as we preferentially detect 
stars moving toward us and are less sensitive to objects 
moving away from us, deeper into the cloud. 

To determine the difference between the radial velocity 
of each protostar and its local molecular gas, we relied on 
measurements from the CO survey of the Milky Way by 
iDame et al.l (|200lD . This survey homogeneously mapped 
all three star forming regions from which the protostars 
studied here were selected, sampled with a 0.125° beam 
(FWHM profile) and a grid spacing (measured in units of 
right ascension and declination) of 0.25° in the Serpens 
and p Ophiuchus molecular clouds, and 0.125° spacing 
in Taurus. The central molecular gas velocity associated 
with each protostar in the sample was measured from the 
profile of the CO line at the nearest sampled position on 
the sky, usually within 1 beam-width of the protostel- 
lar position. CO linewidths were typically 1-2 km sec -1 . 
The central molecular gas velocity was found by bisect- 
ing the portion of the CO line which lies above 25% of 
the peak CO emission at that location. This method 
allows for a robust determination of the central cloud 
velocity which is insensitive to the precise line profile 
of the g as , a desirable quality given the optically thick 
and non-gaussian nature of many of the CO line profiles 
we examined. Typical differences between line center 
velocities given by this method as opposed to standard 
gaussian fits to the CO profile are 0.2 km sec -1 , with 
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FlG. 2. — Histogram display of protostar residual velocities (solid 
line; V stellar - Vco)- The best fit gaussian to this distribution, binned 
in the same manner, is shown as a dashed histogram. The mean and 
a of the gaussian fit arc displayed as well; though the protostars ap- 
pear to have a wider velocity dispersion (cr ~ 2.59 km scc^ 1 ) than 
the standards shown in Figure 1, the rotational broadening and ac- 
cretion veiling present in protostcllar spectral features lead to more 
uncertain velocity determinations. The expected magnitude of these 
effects match well the excess width of the protostcllar velocity disper- 
sion, suggesting we have yet to measure the true velocity dispersion of 
protostcllar populations (see Section 3). 

maximum deviations up to 0.5 km sec -1 . The central 
CO velocity at the nearest sampled position to each pro- 
tostar is shown in the third column of Table 1, and the 
difference between the CO velocity and the protostellar 
velocity is shown in the fourth column. The fifth column 
shows the distance between the protostellar position and 
the center of the CO sampling position, in units of the 
beam profile (FWHM). 

3. THE PROTOSTELLAR RADIAL VELOCITY DISPERSION 

Given the high degree of homogeneity in both the CO 
and NIRSPEC surveys, we are able to examine the radial 
velocity dispersion about the local CO systemic velocity 
of these protostars in bulk. Shown in Figure 2 as a solid 
line is a histogram of the residual protostellar velocities 
about the central velocity of the local CO g as \ a dashed 
line indicates the best fit gaussian to this distribution. A 
comparison of the gaussian fits to the radial velocity dis- 
tributions in Figures 1 and 2 show the protostellar pop- 
ulation has a measured velocity dispersion almost twice 
that measured for our set of standards (standards a — 
1.28 km sec -1 , protostellar a = 2.59 km sec -1 ). 

The excess width of the protostellar velocity dispersion 
above that of the velocity dispersion of the radial veloc- 
ity standards is likely due to factors which make the pro- 
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tostellar velocity determinations less reliable than those 
derived for the standards. The first of these factors is 
the shallower and broader spectral features detected in 
protostellar spectra. The presence of continuum veiling 
and moderate rotational velocities in these protostars re- 
sult in spectral features which are shallower and broader 
than features in standards of similar spectral types. The 
decreased strength of these features result in cross cor- 
relation functions whose central velocity is more uncer- 
tain, demonstrated by examining the standard deviation 
of velocities determined for the same object from differ- 
ent templates; radial velocities determined for standards 
from each of three appropriate templates show a median 
standard deviation of 0.8 km sec -1 , while the median 
standard deviation of the velocities derived for the pro- 
tostellar sample is 1.2 km sec" 1 , 50% larger. Additional 
width may be introduced into the protostellar velocity 
dispersion given the uncertainty in the CO velocity to 
which each protostar is being compared - while our ob- 
servational uncertainty has been examined by comparing 
our derived velocities to known velocities with typical er- 
rors on the order of tens of meters per second, errors in 
the local CO velocity could be as large as 0.2 km sec -1 . 

The combination of weaker protostellar features and 
uncertainty in the local CO velocity implies the obser- 
vational limits of the measured protostellar velocity dis- 
persion are ~ 0.6 km sec -1 larger than of our standards, 
such at a zero velocity population would likely have a 
measured velocity dispersion of ~ 1.9 km sec -1 . This re- 
duces the excess detected width of the protostellar veloc- 
ity dispersion to 0.7 km sec -1 above the level of expected 
observational uncertainties. This excess width could be 
explained by an intrinsic protostellar velocity dispersion 
of 1.8 km sec -1 added in quadrature to the estimated 
1.9 km sec -1 observational uncertainties. It is equally 
plausible, however, that measuring the width of the ve- 
locity dispersion with a gaussian fit to a population with 
31 members will allow additional statistical uncertainties 
due to sampling errors. 

To test the statistical importance of any difference 
between the observed protostellar-CO velocity residuals 
and the observed radial velocity standard residuals, we 
have carried out a Kolmogorov-Smirnov (KS) test on the 
two sets of data. The test is illustrated in Figure 3, 
which compares the cumulative probability functions of 
the protostellar and standard residuals. The KS statistic 
measures the maximum difference between two cumula- 
tive distribution functions and can be used to calculate 
the probability that two distributions were drawn from 
the same parent sample ijPress et al.lll986T) . As seen in 
Figure 3, the two distributions have a nearly 50% chance 
of being derived from the same parent population. 

We therefore find that the dispersion of protostellar ra- 
dial velocities about the local CO gas velocity does not 
differ significantly from the same quantity as calculated 
for a sample of standards with accurately known radial 
velocities. Given our observational errors, this constrains 
the width of the radial velocity dispersion of Class I/flat 
spectrum protostars to below ~ 2.5 km sec -1 . Within 
the limits of our measurements, therefore, we find no evi- 
dence that protostars and their local gas inhabit different 
dynamical environments - that is, at the ~ 2.5 km sec -1 
level, protostars and the local gas appear to be respond- 
ing to the same gravitational potential. Any effect of 
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FlG. 3. — Comparison between the cumulative probability functions 
of the standard and protostellar radial velocity residuals. The proba- 
bility the two functions originated from the same parent distribution, 
as estimated by an application of the KS test, is shown to be 48%. 

stellar interactions on the dynamical state of protostars 
must therefore lie below the ~ 2.5 km sec -1 level, or 
must be dominant only on timescales longer than that of 
the Class I/Flat Spectrum lifetime. 

4. RADIAL VELOCITY OUTLIERS 

Although the radial velocity distribution of the pro- 
tostars is consistent in bulk with that of the surround- 
ing CO gas, we do detect four outliers (04016+2610, 
04489+3042, GV Tau S & CRBR 12) separated from 
the central velocity of the local gas by more than 7.5 
km sec -1 . Given our ±2 km sec -1 radial velocity un- 
certainty, these objects represent 3cr outliers from the 
central velocity of the local CO gas. All four outliers 
possess negative (blueward) velocity shifts with respect 
to the local CO gas and other protostars in the same 
star forming region - the outliers possess anomalous ve- 
locities, and are not outliers due to peculiar local CO 
velocities. The three most significant outliers reside in 
the Taurus- Auriga star forming region. These velocity 
outliers also appear to be bona fide protostars, and not 
older stars along the line of sight to the cluster: all 4 
radial velocity outliers were found by D04 to have r^ > 
0.9, all but CRBR 12 show HI Br 7 emission, and two 
(04016+2610, GV Tau S) show H 2 emission. 

There are 3 plausible explanations for such extreme 
velocity discrepancies in the radial velocities of young 
protostars: a) the protostars have had a close gravita- 
tional encounter with another member of the cluster and 
gained energy enough to be ejected from the cluster, b) 
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TABLE 2 

Maximum Orbital Periods for Possible Binaries 



Target name 


Estimated M 


Maximum Orbital Period 


Binary Separation 


Epoch of Obs. 




(Mq) 


(years) 


(AU) 


(MM-DD-YYYY) 


04016+2610 


1.6 


2.08 


2.40 


11-04-2001 


04489+3042 


0.15 


0.06 


0.11 


11-04-2001 


GV-Tau S 


1.6 


5.55 


4.62 


11-06-2001 


CRBR12 


0.75 


10.45 


5.47 


06-20-2003 



the protostar is a member of a binary with a systemic ve- 
locity equal to that of the local CO gas, but was observed 
at a time when its primary's orbital velocity is directed 
in large part along our observational line of sight, or c) 
the velocities derived for these outliers have large uncer- 
tainties that are poorly characterized by our analysis of 
the velocities derived for stars with well known radial 
velocities. 

The last two scenarios presented above can be eas- 
ily confirmed or rejected with additional observations at 
later epochs. It is plausible that the measurements re- 
ported here are significantly uncertain; three of these out- 
liers have rotational velocities in excess of 25 km sec -1 , 
making cross correlation peaks more difficult to centroid, 
and three of these sources have been observed in the 
optical by White & Hillenbrand ( 2004), who find radial 
velocities for these objects (epoch Dec. 1999) that are 
consistent with other members of the Taurus star form- 
ing region. For these reasons, additional observations to 
confirm the radial velocities reported here are necessary 
before any conclusions can be drawn as to the nature of 
these outliers. However, to guide the interpretation of 
the resulting observations, we briefly examine below the 
possible ramifications of results of followup observations. 

If multiple followup observations show in each epoch 
a protostar with a radial velocity equal to that of the 
local CO gas, it is safe to assume that these objects are 
outliers due to uncertainties in the analysis presented 
here (option c). Additionally, if the radial velocities of 
these outliers are observed to change in a statistically 
significant fashion over the course of the follow up ob- 
servations, it would be likely that their velocity excesses 
result from motions due to binarity (option b). However, 
multiple followup observations showing radial velocities 
consistent with those reported here would not necessar- 
ily be proof that the stars in question are being ejected 
from the star forming region, as it is conceivable that the 
outliers could be members of a binary with a long orbital 
period. 

To constrain the timescales over which a consistent ra- 
dial velocity in one of the outliers would need to ob- 
served before its possible binary status can be ruled out, 
we have constructed a simple model of each outlier as a 
binary system. The binary parameters actually allowed 
for each source by the observed spectra depend on the 
specific T e ff, veiling, radius, and extinction (and spatial 
structure thereof) associated with each source. In gen- 
eral, the observed veiling and lack of significant HI Br 
7 absorption rule out the presence of hot, massive early 
type companions. Accordingly, to provide upper limits 
on the orbital timescales these possible binaries should 
vary on, we will model each system as consisting of equal 



mass components - generating the same radial velocity 
offsets with lower mass companions require smaller or- 
bital separations and shorter orbital timescales. 

Assuming a canonical age for these outliers of 1 Myrs, 
the effective temperatures determined by D05 for the 
sy stem can be t ransla ted into a mass via the models 
of iBaraffe et alJ (^98). Table 2 displays the estimated 
masses derived from the Baraffe models for each of these 
outliers - in our simple model of these outliers as binary 
systems, we adopt the mass estimate given in Table 2 for 
each component of the binary. The two hottest veloc- 
ity outliers, GV Tau S and 04016+2610, have tempera- 
tures ~ 250 K hotter than the hottest stars computed by 
Baraffe et al., requiring a small extrapolation to extend 
the grid to derive their likely masses. 

With mass estimates for the binary system compo- 
nents, we can now calculate the maximum timescale for 
changes in the observed radial velocity if the velocity 
separation from the CO velocity is due to a binary or- 
bit centered on the mean CO gas velocity. To calculate 
this maximum orbital period, we assume the observed 
systems consist of equal mass, single line spectroscopic 
binaries in circular, edge on orbits with orbital velocities 
equal to the velocity separation between the observed ve- 
locity and the local CO gas (shown in column 4 of Table 
1). Under these assumptions, the system's orbital period 
is related to the masses of the stars by: 

P 

com P _ —— v 3 b (1) 

{m obs + m comp ) 2 2ttG 

where m f, s and v b s indicate the estimated mass and de- 
tected residual velocity of the observed component, and 
m comp indicates the upper limit to the mass of any com- 
panion. Our assumption of equal mass binaries implies 
that m b s = m comp . With the observed orbital velocity 
fixed, the assumptions of an equal mass binary in a circu- 
lar, edge on orbit with maximal projection of the orbital 
velocity along the radial velocity direction drive the de- 
rived orbital period to the maximal end of the allowed 
range. Relaxing those assumptions would predict shorter 
orbital timescales and more rapid radial velocity varia- 
tions. Thus, the timescales we derive for the completion 
of a full orbit represent upper limits - changes in radial 
velocity may happen on shorter timescales, but should 
not occur on longer timescales unless the system has a 
primary with T e // greater than that detected by D05. 
The resultant maximum orbital period for each outlier is 
shown in column 3 of Table 2, with corresponding binary 
separations given in column 4. Typical orbital timescales 
for systems under these assumptions are on the order of 
several years, with typical binary separations of 2-5 AU. 
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5. DISCUSSION 

The radial velocity measurements analyzed here repre- 
sent the largest sample of radial velocity measurements 
for Class I/flat spectrum sources yet achieved. These 
measurements allow investigation of a number of dynam- 
ical aspects of low mass star formation. 

As shown in Section 3, the observed protostellar radial 
velocity dispersion (a = 2.59 km sec -1 ) is not signifi- 
cantly broader than would be expected given the obser- 
vational limits of this study. These observations therefore 
place constraints on the width of the protostellar veloc- 
ity dispersion to lie below the level of our observational 
uncertainties, or ~ 2.5 km sec -1 . Velocity dispersions 
of CO gas (as measured using gaussian fits) in these star 
formation regions are typically ~ 1.4 km sec -1 ; the upper 
limit placed on the protostellar velocity dispersion in this 
study is consistent with the local CO gas and protostars 
being dynamically linked during the Class I/flat spec- 
trum stage, with little dynamical evolution from stellar 
interactions. This upper limit does not, however, rule 
out possible dynamical effects which could inflate the 
protostellar velocity dispersion above that of the local 
CO gas, but below the 2.5 km sec -1 level. Future ob- 
servations of protostars producing spectra with higher 
signal-to-noise ratios or more simultaneous wavelength 
coverage than obtained by D05 may allow for more ac- 
curately measured velocities via cross correlation with 
spectral templates, placing tighter constraints on the dy- 
namical state of protostars. Radial velocity uncertainties 
introduced by relatively large rotational velocities, how- 
ever, may prove difficult to overcome. 

Four of the 31 (12.9%) Class I/flat spectrum sources 
appear to have radial velocities more than 3 a away 
from the local cloud velocity. However, given these 
high usinis measured for these sources and their non- 
disc repant radial velocities as de rived from optical spec- 
tra ijWhite fc Hillenbrand! 12004) . followup spectroscopy 
to confirm these radial velocities should be taken before 
these results are overinterpreted. If real, these detections 
may present a lower limit t o the close binary fraction 
of pro tostars. lHaisch et alJ l)2004|) and iDuchene et al.1 
( 2004) find companion star fractions of Class I/flat spec- 
trum objects similar to that of more evolved T Tauri 
stars (~ 20 — 25%) in the star formation regions studied 
here. The photometric nature of their study makes them 
most sensitive to companions with differences in K band 
apparent magnitude less than 4 to 6 magnitudes and at 
orbital distances from 300-2000 AU. A radial velocity 
survey such as this is biased towards finding binaries at 
much closer radii: an equal mass binary system with 0.5 
Mq components in an edge-on, circular orbit with an or- 
bital velocity of 7.5 km sec -1 (the 3 a limit of this study) 
will have a separation of 3.9 AU and a period of 7.8 years. 
Protostars with larger primary masses will have greater 
separations for similar orbital velocities. Thus, the 12.9% 
of stars detected in this work as radial velocity outliers 
could represent a single epoch lower limit to the number 
of protostars with close companions within ~ 10 AU. 

The companion star fraction lower limit presented 
above, however, is uncorrected for incompleteness. As 
part of a multi-epoch radial velocity survey for short 
period (P < 100 days) T Tauri multiple systems, iMelol 
( 2003) calculates incompleteness corrections for binaries 



with a 1.0 Mq primary as a function of period (from 1 to 
1000 days) and secondary mass. Their typical detection 
efficiency is about 60%, though dropping significantly to- 
wards longer periods. Applying a similar correction to 
the number of radial velocity outliers seen here would im- 
ply a true companion star fraction of ~ 22%. This would 
be roughly consistent with the companion star fraction 
measured for T Tauri stars over th e same range of separa- 
tions (see Figure 6 of lMeloll2003L interpolating over the 
gap in data for systems with periods between 100 and 
10000 days ). However, as this study is based on single 
epoch observations, the completeness fraction is likely 
significantly lower than 60%, implying an even higher 
'true' companion star fraction, in excess of that observed 
for T Tauri stars at the same separations. Given the 
good agreement between the T Tauri and protostellar 
companion star fractions at spatially resolvable separa- 
tions, such a difference in multiplicity characteristics at 
close separations would be somewhat surprising. As a 
result, it seems unlikely that all 4 radial velocity outliers 
are short period protostellar binaries. 

The four radial velocity outliers can also be inter- 
preted as single stars with large intrinsic space veloci- 
ties. Such objects have presumably been ejected from 
their initial formation sites by dynamical interactions 
with other cluster memb ers, as suggested for PV Ceph by 
iCoodman fc Arcei l|2004tL The blueward velocity shift of 
all 4 outliers could also be expected as a result of an ejec- 
tion model, as introduced in Section 2.2. However, the 
high velocities reported here (7.5 km sec -1 and larger) 
are unusual even within the context of ejection models, 
as simulations of stellar encounters in a cluster environ- 



ment (|Delgado-Donate et all2003 IrBate et al.ll2UU;-if) pre- 
dict one dimensional velocity dispersions on the order of 
a few kilometers per second, and maximum ejection ve- 
locities of ~ 5 km sec -1 . The detection of spectroscopic 
accretion signatures also suggest the mass reservoir of 
these protostars have not been disrupted recently by dy- 
namical interactions. We therefore view interpretation 
of these high velocity outliers as a population of ejected 
cluster members as unlikely. 

However, our upper limit of ~ 2.5 km sec -1 on the 
one dimensional radial velocity dispersion of protostars 
is not inconsistent with predictions of models of star for- 
mation which include lower velocity dynamical effects. 
As noted earlier, dynamical simulations of clustered star 
formation find 1-D radial velocity dispersions of 1-2 km 
sec -1 , comparable to the level of observational uncertain- 
ties in this study. Additionally, halos of X-ray sources 
have been detected near sites of ongoing star forma- 
tion in ROSAT X-ray surveys. These halos typically 
extend to ~ 5-10 deg rees away from the center of the 
star formation region Feieelsonll 19961: iNeuhaeuser et all 
119071 IWir.hma,nn et a,lJl2000D . The X -ray sources popu- 
lating these halos possess Li absorption features, lead- 
ing some to identify them as a dispersed population of 
Weak lined T Tauri Stars, while others find a lack of 
M type members of this population (which deplete their 
lithium on faster timescales than G & K stars) indicative 
of an olde r population unassociated with ongoing star 
formation (|Briceno et alJll997l iBriceno et al.lll999l) . We 
find that the formation of dispersed halos of young X-ray 
sources could result naturally from the expansion of clus- 
ters whose radial velocity dispersions lie somewhat below 
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the up per limits presented here and bv lHartmann et alJ 
(1986) for T Tauri stars. For instance, at the distance of 
Taurus (140 pc; IMiasf l978) a star moving 1 km sec 1 in 
the plane of the sky can travel 7°. 5 in 18 Myrs, consistent 
with the size and estimated ages of X-ray halos. Thus, 
a definitive analysis of the ability of pre-main sequence 
stars to produce diffuse halos of X-ray sources will require 
stronger constraints on the radial velocity distribution of 
star forming regions than are currently available. 

6. CONCLUSIONS 

Our analysis of the radial velocities of Class I and flat 
spectrum protostars from three nearby star formation 
regions results in the following conclusions: 

• The one dimensional radial velocity dispersion of 
Class I/flat spectrum protostars is measured to be 
<2.5 km sec -1 . 

• The protostellar velocity dispersion upper limit 
presented here is consistent with the measured ve- 
locity dispersions of both CO gas and more evolved 
T Tauri stars in the same star formation regions. 
This is consistent with the dynamics of the CO gas, 
protostars and T Tauri stars being dominated by 
the form of the large scale gravitational potential. 
However, the 2.5 km sec -1 protostellar radial ve- 
locity dispersion upper limit measured here is not 
in itself inconsistent with models of star formation 
in which dynamic stellar interactions play a signif- 
icant role. 

• No large scale bias towards negative radial veloci- 
ties is detected in this sample, as may be expected 
for stars with ejection velocities carrying them to- 
ward Earth and to areas of lower extinction. The 
lower extinctions seen by these members would bias 
magnitude limited samples towards their inclusion, 
thus translating into a stellar median velocity offset 
from that of the local CO gas. 



• We detect 4 of the 31 objects (12.9%) studied here 
as 3 a outliers from the local CO velocity. Further 
observations are needed to confirm these supris- 
ing velocities, which if true, imply either a large 
(< 25%) close (r < 10 AUs) binary fraction, or 
a significant level of dynamical ejections of young 
stars whose accretion reservoirs are nonetheless rel- 
atively undisturbed. The 4 outliers detected here 
do have negative measured radial velocities, con- 
sistent with the blueward bias expected for ejected 
objects. 
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